Submission to Phys. Rev. B 



Magnetic excitations in the spin-5/2 antiferromagnetic trimer 

substance SrMn3P40i4 

Masashi Hase^|*| Masaaki Matsuda^'^, Koji Kaneko^, Naoto Metoki^, Kazuhisa Kakurai^, 
Tao Yang"^, Rihong Cong^, Jianhua Lin^, Kiyoshi Ozawa^, and Hideaki Kitazawa^ 
^National Institute for Materials Science (NIMS), 
1-2-1 Sengen, Tsukuha, Ibaraki 305-0047, Japan 
^ Japan Atomic Energy Agency (JAEA), 
2-4 Shirakata Shirane, Tokai, Naka, Ibaraki 319-1195, Japan 
^Oak Ridge National Laboratory (ORNL), Oak Ridge, TN 37831, USA 
'^College of Chemistry and Molecular Engineering, 
Peking University, Beijing 100871, People's Republic of China 
(Dated: October 3, 2011) 

Abstract 

A quantum-mechanical 1/3 magnetization plateau and magnetic long-range order appear in the 
large-spin (5/2) substance SrMn3P40i4. Magnetization results of SrMn3P40i4 can be explained by 
the spin-5/2 isolated antiferromagnetic linear trimer with the intra-trimer interaction (Ji) value of 
4.0 K. In the present study, to confirm the spin system, we performed inelastic neutron scattering 
(INS) experiments of SrMn3P40i4 powders. We observed plural magnetic excitations. The peak 
positions are 0.46, 0.68, and 1.02 meV. Constant-Q-scan spectra at several Q values (magnitude 
of the scattering vector) indicate that the dispersion is weak. The weak dispersion indicates that 
the excitations are transitions between discrete energy levels. Our INS results are consistent with 
results expected in the trimer model. We evaluated the Ji value as 0.29 meV (3.4 K) without 
considering the other interactions. 
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I. INTRODUCTION 



Quantum-mechanical nature is sometimes apparent even in an ordered state of several 
low- dimensional spin systems formed by small spins. In the triangular antiferromagnet 
CsCuCla with spin-1/2, a small jump was observed in the magnetization curve in the mag- 
netic field parallel to the c axis.|l| This jump was successfully explained as a spin flop 
process caused by ,uaMum-,„echa.ica. effects. H sp„.gap systems with spi„-l/2 such as 
the spin-Peierls system in CuGeOs [3|, |J] and the two-leg ladder system in SrCu203,[5, 6|| 
antiferromagnetic long-range order (AF-LRO) appears when small amounts of impurities 
were doped. |7h11|. Nonetheless excitations originating in the singlet-triplet gap in the pure 



system were observed. 
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It has not been studied adequately whether the quantum-mechanical nature can remain 
in an ordered state of spin systems formed by large spins. Quantum-mechanical nature is 
apparent in the ordered state of the spin-5/2 substance SrMn3P40i4.[l3| The temperature 
T dependence of the magnetic susceptibility indicates occurrence of a magnetic LRO below 
about 2.6 K. A 1/3 magnetization plateau in magnetization curves was observed at both 
1.3 K and 4.2 K (below and above the transition temperature) . jisl. [l^ The plateau was 
observed in a powder sample. Therefore, the plateau can appear irrespective of the applied 
magnetic field direction. As a result, the plateau is a quantum- mechanical magnetization 
plateau generated by an energy gap in the magnetic excitation spectrum. 

It is important to determine the origins of the quantum-mechanical nature (magnetization 
plateau) in SrMn3P40i4. We could explain well the magnetic-field H and T dependences of 
the magnetization using the spin-5/2 isolated linear trimer formed by the AF Ji interaction 
(Ji = 4.0 K) depicted as ellipses in Fig.l.[l^ Therefore, we consider that the magnetization 
plateau originates in discrete energy levels of the AF trimer. If the AF trimer model is 
valid, we can observe magnetic excitations indicating the discrete energy levels even in 
powder samples. Consequently, we performed inelastic neutron scattering experiments of 
SrMn3P40i4 powders. 
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FIG. 1: (Color online) 



5/2 in SrMn3P40i4. 



A schematic drawing of positions of Mn^+ ions (3d^) having localized spin 



13 | Two crystallographically independent Mn^^ sites (Mnl and Mn2) exist. 
Two kinds of short Mn-Mn bonds exist and have Mn-O-Mn paths. The Mn-Mn distances are 3.27 
and 3.34 A at room temperature. The exchange interaction parameters are respectively defined as 
Ji and J2. The dominant AF Ji interactions form the spin trimers indicated by the ellipses. The 
Hamiltonian is expressed as ?^ = Ji{SiS2 + S2S3). The spin trimer can account for the magnetic- 
field and temperature dependences of the magnetization when Ji = 4.0 K. Mn-Mn distances in the 
other bonds are more than 4.89 A. These bonds have no Mn-O-Mn paths. 



II. METHODS OF EXPERIMENTS 



We synthesized single crystals of SrMn3P40i4 under hydrothermal conditions at 473 
K. jisjl Each crystal was small. We used pulverized crystals for inelastic neutron scattering 
(INS) measurements. 

We carried out INS measurements on the cold neutron triple-axis spectrometer LTAS 
installed at JRR-3M in JAEA. The final neutron energy was fixed at 2.6 meV. Higher-order 
beam contamination was effectively eliminated using a cooled Be filter before the sample. 
The horizontal collimator sequence was guide-80'-Be-sample-120'-open. This setup yields an 
energy resolution of 0.1 meV (full width at half maximum, FWHM) at an energy transfer 
oj = meV. The resolution was determined from incoherent scattering of the sample. The 
powder sample of about 9 g was mounted in a ^He closed cycle refrigerator. 



III. RESULTS AND DISCUSSION 



We performed all the INS measurements above the transition temperature. Circles in 
Fig. 2 show the u dependence of the INS intensity (constant-Q scan spectra) around 5 K. 




<2;(meV) 



FIG. 2: (Color online) The inelastic neutron scattering intensity vs. energy transfer (w) [constant- 
Q (magnitude of the scattering vector) scan spectra] of SrMn3P40i4 around 5 K (circles). The 
solid line represents the sum of two Gaussians and one Lorentzian (plus constant backgrounds). 
The dashed line indicates each Gaussian or Lorentzian. The values of the parameters are given in 
Table I. 

The value of Q is the magnitude of the scattering vector. Excitations are apparent between 
0.5 and 1.5 meV. Two kinds of excitations with different peak positions seem to overlap each 
other. The spectra are almost independent of Q except for difference in intensities. The 
weak Q dependence indicates that the excitations are transitions between discrete energy 
levels. We consider that the intensities in the vicinity of meV cannot be explained only 
by incoherent scattering because of T dependence of spectra as shown later. Low-energy 
excitations exist. 
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We compared each spectrum above 0.2 meV in Fig. 2 with a sum of two Gaussians and 
one Lorentzian (plus constant backgrounds). 

m - + exp[-i^^^] + 7bg. (1) 

Here the sum is from z = 1 to 2. The two Gaussians correspond to the excitations between 
0.5 and 1.5 meV. The Lorentzian corresponds to the excitations in the vicinity of meV. 
Each sum of two Gaussians and one Lorentzian (sohd hne) reproduces well the corresponding 
spectrum in Fig. 2. Obtained values of the fitting parameters are shown in Table L The 
peak position in the spectrum at Q = 1.0 is 0.68 meV or 1.02 mcV. The peak width 
(FWHM) is 0.28 meV and 0.55 meV for the 0.68 meV and 1.02 meV excitation, respectively. 
These widths are larger than the energy resolution of 0.1 meV at = meV, indicating 
existence of weak dispersion caused by inter-trimer interactions. 

Circles in Fig. 3 show constant-Q scan spectra at Q = 1.0 A~^. The 0.68 meV and 1.02 
meV excitations are also seen at 11.2 K. Intensities around 0.5 meV are larger at 11.2 K 
than at 5.0 K, suggesting appearance of another transition. Therefore, we compared the 
spectrum above 0.2 meV at 11.2 K with a sum of three Gaussians and one Lorentzian (plus 
constant backgrounds) given in Eq. (1) with i = 1 to 3. To reduce variable parameters, we 
assumed that the peak position (0 meV) of the Lorentzian, and the peak positions (0.68 and 

I. 02 meV) and widths of the two Gaussians were constant. We used the values obtained 
at 5.0 K. This assumption is reasonable for transitions between discrete energy levels. The 
sum of the three Gaussians and one Lorentzian (solid line) reproduces well the spectrum at 

II. 2 K. The peak position of the third Gaussian is 0.46 meV. The peak width (FWHM) is 
0.25 meV and is larger than the energy resolution of 0.1 meV at a; = mcV. 

The spectra at 15.6 K and 20.3 K shown respectively in Figs. 3(b) and (c) resemble the 
spectrum at 11.2 K. Therefore, we compared the spectrum above 0.2 meV at 15.6 K or 20.3 
K with a sum of three Gaussians and one Lorentzian (plus constant backgrounds) given in 
Eq. (1) with i = 1 to 3. In the fitting, we assumed that the peak position of the Lorentzian, 
and the peak positions and widths of the three Gaussians were constant. The sum (solid 
line) reproduces well each experimental spectrum. Spectra are featureless above 30 K. We 
did not compare the spectra with calculated curves. The integrated intensity between 0.5 
and 1.5 meV decreases slightly on heating. The Bose factor proportional to phonon intensity 
at 0.68 meV and 20.3 K, on the other hand, is about 9 times as large as that at 0.68 meV 
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TABLE I: Values of the integrated intensity (Ij) and full width at half maximum (FWHM) of 
Lorentzian or Gaussian obtained from the fitting of Eq. (1) to experimental constant-Q spectra of 
SrMn3P40i4. FWHM is given as 2ao for Lorentzian or 2\/ln 2aj for Gaussian. In the fitting to the 
spectra at Q = 1.0 A~^, the value of FWHM of each Gaussian was obtained at 5.0 or 11.2 K and 
was fixed in the fitting at higher T. The values in parentheses indicate errors. 
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and 5.0 K. Therefore, contribution of phonon is small enough and magnetic excitations are 
dominant between 0.5 and 1.5 meV. 

We examined whether the spin-5/2 AF trimer model with Ji = 4.0 K can account for the 
observed excitations. Figures 4(a) depicts a schematic drawing of low-lying energy levels, [l^ 



The following selection rules of transitions are derived theoretically. 15(] 



AS = 0,±1, AM = 0,±1, and A^i3 = 0,±1. (2) 

Si{i = 1,2, 3) is the spin operator in the trimer. S and Srs is defined as Si + 5*2 + S3 and 
Si + S3, respectively. M is the z component of S. Arrows in Fig. 4(a) indicate allowed 
transitions from the ground state (GS), first-excited state (lES) or second-excited states 
(2ES). In our experimental setup, we can observe transitions with an energy difference Ae 
up to 8 when Ji = 4.0 K. Here, e is defined as E/Ji [E: eigen energy). 
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FIG. 3: (Color online) Constant-Q scan spectra of SrMn3P40i4 at Q = 1.0 (circles). The 
solid line represents the sum of three Gaussians and one Lorentzian (plus constant backgrounds). 
The dashed line indicates each Gaussian or Lorentzian. The values of the parameters are given in 
Table I. 



Figure 4(b) depicts T dependence of calculated occupation ratio of the five low-lying 
energy levels. An inelastic neutron scattering intensity strongly depends on the occupation 
ratio. From Fig. 4(b), we know that excitations from GS are dominant around 5 K. We 
considered that the 0.68 and 1.02 meV excitations correspond to transitions from GS to lES 
and 2ES, respectively, indicated by black arrows. The respective energy differences are 2.5 Ji 
and 3.5Ji. The value of Ji is evaluated as 0.27 meV (3.2 K) or 0.29 meV (3.4 K). These 
values are slightly smaller than the value determined in the magnetization results ( Ji = 4.0 
K). Excitations from lES or 2ES are also expected at 11.2 K. We considered that the 0.46 



7 



8ES, *=-6.5,S=7/2, 5,, = 4 — 
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FIG. 4: (Color online) (a) A schematic drawing of lo w-ly ing energy levels in the spin-5/2 AF linear 



l4 \ The parameters e = E/Ji and S indicate 



trimer [ground state (GS) and excited states (ESs)]. 
the eigen energy and the total spin, respectively. S'13 is defined in the text. To distinguish two 
degenerate eigen states with e = —10, we name the two states 3ES and 4ES. To distinguish two 
degenerate eigen states with e = —7, we name the two states 6ES and 7ES. The arrows indicate 
allowed transitions from GS, lES or 2ES with an energy difference Ae up to 8. We observed the 
three black-arrow transitions. The gray-arrow transitions may exist. We could not detect the 
white-arrow transition, (b) The temperature dependence of calculated occupation ratio of the five 
low-lying energy levels when Ji = 4.0 K. 



meV excitation corresponds to the transition from 2ES to 4ES indicated by a black arrow. 
The energy difference is 1.5Ji. The value of Ji is evaluated as 0.30 meV (3.5 K) and is close 
to the values evaluated from the other two transitions. 

We examined whether we observed all the allowed transitions that are possible in our 
experimental setup. As was described, we observed the three black-arrow transitions. The 
gray- arrow transitions may exist. However, we could not prove the existence of the gray- 
arrow transitions. Energy differences of some gray- arrow transitions (2.5Ji and 4Ji) are 
the same as or close to the energy differences of black- arrow transitions (2.5Ji and 3.5Ji). 
Therefore, we could not extract contribution of the gray-arrow transitions from the experi- 
mental results. Energy differences of the other gray-arrow transitions are 4.5Ji = 1.3 meV 
and 5Ji = 1.4 meV when Ji = 0.29 meV. Small INS intensities are seen around these ener- 
gies in Figs. 2 and 3. In our analyses, the small intensities correspond to the tail of the 1.02 
meV excitation. However, the peak width is larger in the 1.02 meV excitation (FWHM = 
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FIG. 5: (Color online) The temperature dependence of the integrated intensity of the 0.46, 0.68, 
and 1.02 meV excitations in SrMn3P40i4 ai Q = 1.0 A~^. 

0.55 meV) than in the 0.46 meV excitation (FWHM = 0.25 meV) or the 0.68 meV excitation 
(FWHM = 0.28 meV). The gray-arrow transitions may exist in the tail. We could not detect 
the white-arrow transition. We do not have theoretical INS intensities. Therefore, we could 
not determine the reason why we could not detect the white-arrow transition. The INS 
intensity of the white-arrow transition may be very small. Transitions from 3ES or higher 
excited states must exist at 11.2 K and higher T. However, we could not prove existence of 
these transitions because of the same reason for the gray-arrow and white-arrow transitions 
in Fig. 4(a) 

Figure 5 shows the T dependence of the integrated intensity of the 0.46, 0.68, and 1.02 
meV excitations. The integrated intensity of the 0.46 meV excitation is nearly independent 
of T. The integrated intensity of the 0.68 or 1.02 meV excitation gradually decreases with 
increasing T. As is shown in Fig. 4(a), several transitions are expected to exist. Therefore, 
the T dependence of the integrated intensity in Fig. 5 cannot be compared directly with 
the occupations ratios in Fig. 4(b). 

The Q dependence of the INS intensity in the AF trimer is given in the following 



formula. 
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17| 



HQ) = AJ{Qf[l - sin(3.27g)/(3.27g)] + A2f{QY[l - sin(6.54Q)/(6.54g)]. (3) 



The values 3.27 and 6.54 indicates the Mnl-Mn2 and Mn2-Mn2 length in the AF trimer, 
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FIG. 6: (Color online) (a) The Q dependence of the calculated INS intensity in the spin-5/2 AF 
linear trimer. (b)-(d) The inelastic neutron scattering intensity vs. Q (constant-w scan spectra) 
of SrMn3P40i4 (circles). The line indicates the first term in Eq. (3). 

respectively. The function f{Q) is the magnetic form factor of Mn^"*" ions. The coefficients 
Ai and A2 depend on two eigen states between which the transition occurs. The coefficients 
are not derived theoretically. Figure 6(a) represents the two terms in Eq. (3). Circles in Figs. 
6(b) - (d) show constant-w scan spectra of SrMn3P40i4. The lines indicate the first term in 
Eq. (3) plus constant backgrounds. The INS intensity in the vicinity of a; = 2 meV is small 
and almost independent of co, Q, and T. Therefore, we used the intensity at w = 2 meV 
for the value of constant backgrounds. The each line is consistent with the corresponding 
constant-w scan spectrum. If Ai is much larger than A2 in the observed transitions, this 
consistency indicates that the AF trimer model can explain the experimental I{Q). 

We comment on inter-trimer interactions. The dispersion relation of magnetic exci- 
tations was calculated in spin dimers with weak inter-dimer interactions using random 
phase approximation, {lol A similar dispersion relation was inferred in interacting spin 
tetramers. |20[| According to the results, we speculate that the following dispersion relation 



may be applicable to spin trimers with weak inter-trimer interactions. 

uJ^=(h,k,i) = VA2 + «AJ(q)i?(T). (4) 
Here, A is an energy difference between ground and excited states, a is a coefficient derived 
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from transition matrix elements. The value of a is 2 for the spin-1/2 dimer 21 1 or 5 for 

n 

the spin-3/2 dimer. [19| J(q) is a Fourier transform of inter-trimer interactions. -R(T) is 
a difference in thermal populations of ground and excited states. We consider that the 
dominant inter-trimer interaction is the J2 interaction. J(q) is expressed approximately as 
2J2Cos(2TTk). We assume that excitation energies are the same at and 5.0 K. The excitation 
energy at K at the bottom of the dispersion ut,, where the INS intensity is the strongest, 
is expressed as follows. 

Wb = - 2«A J2. (5) 

We assume that only the J2 interaction is the origin of the difference between the expected 
excitation energy 2.5 Ji = 10 K and experimental excitation energy 0.68 meV = 7.9 K in the 
transition between GS and lES. Using A = 10 K and Wb = 7.9 K, we obtained aJ2 = 1.9 
K. If a is large, a J2 value can be small enough in comparison with the Ji value. 

We observed magnetic excitations that are consistent with excitations expected in the 
spin-5/2 AF trimer. Therefore, the discrete energy levels of the AF trimer are the origins of 
the quantum-mechanical nature (magnetization plateau) in SrMn3P40i4. In the strict sense, 
the energy difference between GS (S=5/2) and 2ES (S=7/2) generates the magnetization 
plateau. The magnetization plateau appears even in the ordered state. Q, 1^ The property 
of the cluster (trimer in this case) can remain in the ordered state. The total spin of the 
ground state of the AF trimer is finite (5/2). Therefore, we consider that the magnetic LRO 
is stabilized by the Ji and weak three-dimensional inter-trimer interactions. Several cluster 
substances can maintain their cluster properties in their ordered states. For example, the 
1/2 quantum- mechanical magnetizati on p lateau is generated by discrete energy levels of a 



spin-1/2 tetramer in Cu2CdB206. 22|, |23| The plateau remains in the ordered state. With 



the aid of other researchers, some of the present authors determined the magnetic structure 



below the transition temperature of Tn = 2.2(1) K using neutron powder diffraction data. 24 1 
The magnetic structure has a long-range period. We are now considering the origin of the 
long-range period. We will report details of the magnetic structure in a subsequent paper. 

We comment on the INS intensity in the vicinity of meV. Figure 7 shows constant-Q 
scan spectra aX Q = 1.0 A~^ below 0.5 meV. As was described, the increase of the intensity 
around 0.4 meV is caused by the 0.46 meV excitation. The intensity between 0.1 and 
0.2 meV decreases with increasing T. This temperature dependence cannot be explained by 
incoherent scattering or phonon. In addition, we observed diffuse scattering between 26 = 15 
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FIG. 7: (Color online) Constant-Q scan spectra of SrMn3P40i4 at Q = 1.0 A ^. 



and 40° in neutron powder diffraction patterns (wavelength A = 2.458 A). 2J] This 26 range 
corresponds to Q = 0.7 to 1.7 A^^ in the present INS experiments. The shape of the diffuse 
scattering resembles one-dimensional or two-dimensional Bragg scattering with a cutoff at 
low Q and long tail at large Q. The integrated intensity of the diffuse scattering shows 
a maximum in the vicinity of Tn. Several magnetic reflections appear below Tn between 
26 = 15 and 40°. Therefore, the origin of the diffuse scattering is magnetic. Consequently, 
magnetic excitations exist in the vicinity of meV. The magnetic excitations cannot be 
explained by transitions between energy levels in the trimer. Spin fluctuation in the ground 
state generates the magnetic excitations. Therefore, we used the Lorentzian with the 
meV peak in the fitting of the constant-Q scan spectra in Figs. 2 and 3. In future studies, 
we will perform INS measurements of SrMn3P40i4 in the ordered state. Anisotropy of the 
Mn^"*" spins is small. [iJI Therefore, a gap of spin-wave excitations is small. We may observe 
spin-wave excitations in the vicinity of meV in addition to the trimer excitations. 



IV. SUMMARY 



In order to confirm the spin system, we performed inelastic neutron scattering (INS) 
experiments of powders of the spin-5/2 antiferromagnetic trimer substance SrMn3P40i4. We 
observed plural magnetic excitations. The peak positions are 0.46, 0.68, and 1.02 meV. The 
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weak Q dependence of constant-Q-scan spectra indicates that the excitations are transitions 
between discrete energy levels. The experimental results are consistent with results expected 
in the trimer model with the intra-trimer interaction value of 0.29 meV (3.4 K) without 
considering the other interactions. 
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